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Abstract 

Partial  oxidation  reforming  of  methane  using  a  ceramic  mixed  ionic  (oxygen)/electronic  conducting  dense  ceramic  membrane  within  a 
solid  oxide  fuel  cell  (SOFC)  stack  is  potentially  attractive  as  it  allows  direct  operation  on  methane.  We  have  analysed  this  process  using 
a  simple  model  of  a  counter-flow  fuel  system  in  which  steam  and  carbon  dioxide  from  the  anode  channel  of  the  SOFC  provide  oxygen 
which  passes  through  the  mixed  conducting  membrane  to  partially  oxidise  methane  on  the  other  side  of  the  membrane.  The  model  shows 
that  the  concept  is  feasible  and  that  the  efficiency  is  the  same  as  direct  methane  operation.  The  model  also  gives  target  values  for  the 
properties  of  the  membrane.  Tracer  oxygen  diffusion  and  exchange  have  been  measured  in  the  family  of  oxides  Lai-^Sr^Feo.sCro^Os-g 
under  reducing  conditions.  The  data  show  that  these  materials  have  the  mass  transport  characteristics  and  stability  required  for  this 
concept. 
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1.  Introduction 

Ideally,  fuel  cells  should  be  able  to  run  directly  on  the  raw 
fuel  of  choice,  but  this  is  currently  not  feasible  for  SOFCs  fu¬ 
elled  by  natural  gas  (mainly  methane).  There  has  been  some 
success  with  direct  methane  operation  using  ceria-based 
anodes  [1,2],  but  there  are  still  significant  problems  with 
the  high  anode  over-voltage  and  carbon  deposition.  Conse¬ 
quently,  it  is  usual  to  reform  the  methane  to  H2  and  CO  be¬ 
fore  electrochemical  oxidation  at  the  anode.  Internal  steam 
reforming  on  the  anodes  is  the  most  desirable,  but  the  amount 
of  steam  required  to  suppress  carbon  deposition  is  much 
greater  than  the  1:1  steam:methane  ratio  required  for  the 
reforming  reaction.  The  extra  steam  absorbs  energy  (reduc¬ 
ing  efficiency),  dilutes  the  anode  reactions  (reducing  power 
density)  and  increases  balance-of-plant  complexity  (rais¬ 
ing  steam  and/or  recycling  it).  Furthermore,  the  rapid  en¬ 
dothermic  reforming  reaction  on  the  first  few  anodes  of  the 
stack  gives  a  large  cooling  effect  which  reduces  output  and 
can  cause  mechanical  degradation,  particularly  with  ceramic 
structures  [3].  External  steam  reforming  or  in-stack  reform¬ 
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ing  are  more  controllable,  but  again  require  excess  steam 
and  lead  to  a  loss  in  efficiency  unless  there  is  very  good 
transfer  of  heat  from  the  cell  reaction  to  the  reforming  re¬ 
action.  Partial  oxidation  reforming  requires  no  added  steam, 
but  reduces  efficiency  considerably. 

Several  suggestions  have  been  made  to  overcome  these 
problems.  For  example,  a  concept  of  gradual  internal  re¬ 
forming  has  been  proposed  to  avoid  the  cooling  of  the  first 
anodes  [4].  However,  this  is  still  potentially  susceptible  to 
carbon  deposition.  As  an  alternative,  Steele  [5]  suggested 
using  the  steam  and  CO2  created  by  the  SOFC  cell  reaction 
to  partially  oxidise  the  methane  feed  by  transfer  of  oxygen 
through  a  dense  oxygen-permeable  ceramic  membrane.  The 
transfer  would  simultaneously  boost  the  H2  and  CO  concen¬ 
trations  in  the  anode  stream  and  thereby  increase  the  out¬ 
put  of  the  cells.  Membrane  reforming  of  methane  by  par¬ 
tial  oxidation  of  oxygen  from  air  is  an  area  of  great  interest 
in  its  own  right  and  oxides  with  suitable  properties  are  be¬ 
ing  developed  [6]  that  could  be  adapted  for  application  in 
SOFCs.  The  work  described  here  has  two  objectives.  The 
first  is  to  show  by  a  modelling  approach  that  membrane 
reforming  of  methane  in  a  SOFC  is  feasible  in  principle. 
The  second  is  to  identify  potential  membrane  materials  that 
have  the  required  thermo-chemical  stability  and  transport 
properties. 
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Nomenclature 

As 

area  of  a  reaction  site 

CO 

concentration  of  oxygen  in  solid  phase 

C,  Ca,i.  Q,i 

molar  reformate  concentration  in  ith 
cell  of  anode  or  reformer  channel 

CT 

total  volumetric  molar  gas  concentration 

D* 

oxygen  tracer  diffusion  coefficient 

F 

Faraday’s  constant 

Fa 

gas  velocity  in  anode 
channel 

Fr,i 

gas  velocity  in  ith  cell  of  reformer 
channel 

h 

height  of  reformer  and  anode  channels 

current  densities  in  fuel  cells  and 
reformer  membrane 

jo 

oxygen  flux  density 

k\,  k2 

interfacial  exchange  coefficients 
(rate  constants)  at  membrane  surfaces 

Ka,Kr 

equilibrium  constants  for  gas  in  anode 
and  reformer  channels 

l 

length  of  each  fuel  cell 

L 

membrane  thickness 

Na 

Avagadro’s  number 

No 

number  of  moles  of  oxygen 

Nr 

turnover  frequency  for  reformer  reaction 
with  methane 

Pa,Pr 

oxygen  activities  in  anode  and  reformer 
channels 

-Po2,l>  Po2,2 

oxygen  activities  at  membrane  surfaces 

rc,  rm 

area-specific  resistances  of  fuel  cells 
and  reformer  membrane 

R 

gas  constant 

t 

time 

T 

temperature 

Vmo 

volume  of  oxide  per  gram  atom 
of  oxygen 

Vn 

Nernst  potential 

w 

width  of  each  fuel  cell 

Greek  letters 

y 

thermodynamic  factor  for  oxygen  in 
solid  phase 

MO 

chemical  potential  of  oxygen  atoms 

2.  Transport  through  a  mixed  conducting  ceramic 
membrane 

The  overall  rate  of  permeation  of  oxygen  through  the 
membrane  is  controlled  by  reaction  at  the  two  membrane  sur¬ 
faces  and  the  transport  of  oxygen  ions  and  electrons  through 
the  membrane.  The  membrane  material  must  be  a  mixed 
ionic  and  electronic  conductor  with  a  high  oxygen  ion  con¬ 
ductivity.  For  materials  in  which  the  electronic  conductivity 


is  much  greater  than  the  ionic  conductivity  the  steady  state 
oxygen  flux  density  (oxygen  atoms  per  unit  area  per  unit 
time)  through  a  membrane  of  thickness  L,  is  given  by  the 
expression  [5] 

.  _  _ InUW/JW)  (1-. 

JO  -  2Vmo  (L/(D*))  +  (lAa)  +  (lAc) 


where 


fD*  d(ln  Pq2) 

/acd(lnP02) 


and  Na  is  Avagadro’s  number.  hd  and  kc  are  the  surface 
reaction  rate  constants  (ms-1)  for  the  anodic  and  cathodic 
interfacial  reactions  of  oxygen  at  the  membrane  surfaces 
defined  by  the  following  equation  [7,8]: 


j  (O  atoms  equivalent) 


Na  A/xo 
Vmo  RT 


In  these  expressions  the  subscripts  1  and  2  denote  the  gas 
phases  on  the  two  sides  of  the  membrane  and  a  and  c  are 
the  anodic  and  cathodic  reaction  surfaces  of  the  solid  mem¬ 
brane.  The  activities  of  molecular  oxygen  at  these  interfaces 
are  not  equal  to  those  of  the  corresponding  gas  phases  be¬ 
cause  of  the  chemical  driving  force  of  the  interface  reac¬ 
tions.  This  is  accounted  for  by  the  drop  in  atomic  oxygen 
chemical  potential,  A/xo,  in  going  from  the  gas  to  the  solid. 
(Note  that  A/xo  =  (l/2)A/xo2  for  molecular  oxygen.)  Vmo 
is  the  volume  of  the  oxide  per  gram  atom  of  oxygen  (and  is 
approximately  12.1cm3  for  the  doped  LaFe03  oxides  dis¬ 
cussed  later). 

The  membrane  can  also  be  regarded  as  an  equivalent  elec¬ 
trical  circuit  in  which  the  electro-motive  force  is  the  Nernst 
potential,  Vn»  of  the  two  gaseous  oxygen  activities 


and  the  current  is  carried  by  the  oxygen  ions.  The  equivalent 
area- specific  resistance  of  the  membrane  is  then  given  by 


r  m 


1 


VmoRT 
4  F2 


3.  Membrane  reforming  in  the  SOFC 

In  order  to  evaluate  the  concept  we  have  constructed  a 
simple  SOFC  model  with  an  integrated  ceramic  membrane 
reformer.  The  model  uses  simplified  descriptions  of  the  elec¬ 
trochemical  and  transport  phenomena.  It  is  recognised  that 
in  reality  these  processes  are  much  more  complicated,  but 
nevertheless  they  serve  to  illustrate  the  coupling  between  the 
SOFC  cells  and  the  membrane  reformer.  Fig.  1  shows  the 
model  used.  It  consists  of  20  planar  cells  connected  in  series 
such  that  the  current  is  the  same  in  each  cell  and  the  cell 
potential  is  different  for  each  cell.  This  configuration  is  sim¬ 
ilar  to  the  “segmented  in  series”  design  [9,10].  (Although 


A.  Atkinson,  T.  Ramos /  Journal  of  Power  Sources  130  (2004)  129-135 


131 


Exhaust 

Methane 


<=> 


we  consider  a  specific  SOFC  configuration  in  the  model, 
this  is  only  to  illustrate  the  membrane  reforming  concept. 
The  broad  conclusions  should  be  applicable  to  other  geome¬ 
tries.)  Methane  enters  the  reformer  channel  and  is  reformed 
by  partial  oxidation  by  oxygen  diffusing  through  the  mem¬ 
brane  from  the  anode  channel.  The  partial  reactions  at  the 
two  surfaces  of  the  membrane  are 


H20  +  2e“  -*  O2"  +  H2  and 

C02  +  2e~  O2"  +  CO  (6) 

in  the  anode  channel  and 

CH4  +  O2-  -*■  CO  +  2H2  +  2e~  (7) 


in  the  reformer  channel. 

At  the  end  of  the  reformer  channel  the  reformate  stream 
is  reversed  into  the  anode  channel  of  the  fuel  cells.  The  flow 
of  oxygen  through  each  fuel  cell  into  its  anode  channel  is 
given  by  the  cell  current  density  ic  and  the  cell  potential  of 
each  cell  by 

Tceii  —  Vn,c  —  icr  c  (8) 


where  rc  is  the  area-specific  resistance  (ASR)  of  the  fuel 
cells.  This  assumes  that  the  polarisation  of  the  electrodes  is 
linear. 

The  flow  of  oxygen  across  the  reformer  membrane  is 
driven  by  the  difference  in  oxygen  activities  between  the 
anode  compartment  and  the  reformer  compartment  of  each 
cell.  The  equivalent  electrical  current  is  given  by 


FN,m 

bn  — 

r  m 


where  VN,m  is  the  Nernst  potential  across  the  membrane  and 
the  ASR  of  the  membrane  is  given  by  Eq.  (5). 

The  mass  balance  equations  for  each  cell  include  the  ma¬ 
terial  entering  and  leaving  the  cell  by  advective  flow  and  the 
oxygen  entering  and  leaving  across  the  fuel  cells  and  the 
membrane.  Each  cell  is  assumed  to  have  no  gradients  other 
than  through  its  thickness  and  therefore  a  finite  element  ap¬ 
proach  is  used.  The  moles  of  oxygen  ions  crossing  each  fuel 
cell  per  unit  time  is 


dNo\  wl  . 

~3t~)  c  =  2 Fh 


(10) 


where  w  is  the  width  of  a  cell  and  /  its  length.  The  corre¬ 
sponding  number  crossing  the  MIEC  membrane  is 


wl  1  RT 


IF  rm  4 F 


In 


a 


(ID 


where  Ra  and  Pv  are  the  oxygen  activities  in  the  anode  and 
reformer  sections  of  the  cell. 

Thermodynamic  equilibrium  is  assumed  locally  in  each 
cell.  Gas  equilibrium  calculations  [11]  indicate  that  the  ra¬ 
tio  of  H2  to  CO  should  be  approximately  2  throughout  both 
the  reformer  and  anode  channels  (when  C  deposition  is  sup¬ 
pressed).  Therefore,  in  order  to  simplify  the  model,  it  is  per¬ 
missible  to  treat  the  reformate  mixture  as  though  it  were  all 
H2,  but  at  an  effective  reformate  concentration  equal  to  1.5 
times  the  true  H2  concentration.  Thus  the  mass  balance  for 
the  reformate  concentration  (C  moles  per  unit  volume)  in 
the  anode  channel  of  the  ith  cell  gives 


dCaJ  _  _  r  \  _  *c  1  ^ 

At  ~  l  (Ca,l’+1  aj)  2 Fh  +  8 F2hrm 


In 


a 


(9) 


(12) 
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In  this  equation  Ea  is  the  linear  flow  velocity  in  the  anode 
channel  and  h  is  the  height  of  the  channel.  The  flow  veloc¬ 
ity  is  constant  in  the  anode  channel  because  both  the  fuel 
cell  anode  reaction  and  the  membrane  reaction  conserve  the 
number  of  gas  molecules.  The  first  term  in  Eq.  (12)  is  the  ad- 
vective  contribution  (bearing  in  mind  the  reversed  gas  flow 
means  that  it  is  flowing  from  high  to  low  cell  indices),  the 
second  term  is  consumption  of  reformate  by  the  fuel  cell 
reaction  and  the  third  term  is  increase  in  reformate  concen¬ 
tration  by  loss  of  oxygen  through  the  MIEC  membrane.  The 
corresponding  mass  balance  for  the  reformer  channel  is 


d Crj 

d  t 


l  (Fr,i— 1  Cr,z  —  1  + 


3  RT 


8  F2hrm 


In 


(13) 


The  first  term  is  again  the  advective  contribution,  but  this  is 
slightly  more  complicated  because  the  flow  velocity  in  the 
reformer  channel  is  not  constant  because  the  reforming  reac¬ 
tion  increases  the  number  of  gas  molecules.  The  second  term 
is  due  to  oxygen  permeating  the  MIEC  membrane  and  the 
factor  of  3  is  because  each  oxygen  atom  passing  through  the 
membrane  produces  three  molecules  of  reformate  (Eq.  (7)). 
The  increase  in  reformate  flow  velocity  at  each  cell  is  given 
by 


(14) 


where  P  is  the  total  pressure.  The  anode  flow  velocity  is 
equal  to  the  reformer  flow  velocity  in  the  last  cell  (/  =  20) 
because  the  output  of  the  reformer  is  delivered  into  the  anode 
channel  at  this  cell. 

Gas  phase  equilibrium  in  the  anode  channels  is  determined 
by  the  reactions: 

H2  +  i02  **  H20  (15) 


and 

C0+±02**C02  (16) 


Since  these  have  similar  free  energy  changes  (— 182,944  kJ 
for  Eq.  (15)  and  — 180,424  kJ  for  Eq.  (16)  [11])  then  we  can 
approximate  them  as  being  equal  (this  is  why  the  reformate 
retains  stoichiometric  composition).  Thus  the  equilibrium 
oxygen  activity  in  the  anode  channel  can  be  described  by 
the  following  equation: 


CT  -  C 
K*C 


(17) 


where  Ct  is  the  total  number  of  gas  moles  per  unit  volume 
at  a  pressure  of  1  atm  (10.25  m~3)  and  K&  is  the  equilibrium 
constant  for  Eq.  (15)  (=  1.4  x  108  at  900  °C  [11]). 

Gas  phase  equilibrium  in  the  reformer  channel  is  con¬ 
trolled  by  the  reaction 

CH4  +  j02  CO  +  2H2  (18) 


Assuming  stoichiometric  conversion  then  the  concentration 
of  CO  is  0/3  and  that  of  H2  is  20/3.  Equilibrium  of  Eq.  (16) 
then  gives  the  oxygen  activity  in  the  reformer  channel  as 


4  c3 

21KxC\{Cj 


(19) 


where  Kr  is  the  equilibrium  constant  for  Eq.  (16)  (=1.93  x 
1011  at  900  °C  [11]). 

The  initial  condition  employed  to  give  a  stable  start  was 
that  all  the  reformer  and  anode  channels  were  initially  filled 
with  a  mixture  of  reformate  with  10%  CH4  and  the  CH4 
was  not  allowed  to  react  in  the  fuel  cells.  For  t  >  0  the  input 
to  the  reformer  channel  was  99.9%  CH4  and  0.1%  refor¬ 
mate  at  constant  flow  rate.  The  cathodes  were  assumed  to 
be  in  equilibrium  with  air  at  atmospheric  pressure.  Eqs.  (12) 
and  (13)  were  solved  numerically  by  finite  difference  time 
step  (typically  one-tenth  of  the  residence  time  per  cell)  us¬ 
ing  Maple  software  (Waterloo  Maple  Inc.).  The  simulations 
were  carried  out  for  a  constant  temperature  of  900  °C,  pres¬ 
sure  of  1  atm,  cell  ASR  of  0.5  £2  cm2  and  fuel  utilisation  of 
approximately  85%.  Each  cell  was  assumed  to  have  width 
of  100  mm  and  length  in  the  flow  direction  of  20  mm.  The 
heights  of  the  reformer  and  anode  channels  were  both  taken 
to  be  5  mm.  The  simulations  were  run  until  they  reached 
steady  state.  In  order  to  provide  a  comparison,  a  simulation 
was  also  carried  out  with  stoichiometric  external  steam  re¬ 
forming  (and  no  membrane  reforming). 

The  input  flow  velocity  and  the  MIEC  membrane  ASR 
were  then  adjusted  to  achieve  close  to  full  methane  conver¬ 
sion  at  the  end  of  the  reformer  channel  and  85%  fuel  utilisa¬ 
tion.  This  is  a  limiting  case  since  if  conversion  is  incomplete, 
unreformed  methane  will  enter  the  fuel  cell  anodes,  leading 
to  a  risk  of  carbon  deposition  on  the  anodes.  The  reformate 
profiles  in  the  reformer  and  anode  channels  are  shown  in 
Fig.  2.  The  reformate  fraction  is  defined  as  the  sum  of  the 
mole  fractions  of  H2  and  CO  in  the  gas  stream.  The  model 
demonstrates  that  a  stable  steady  state  exists  to  achieve  ap¬ 
proximately  complete  membrane  reforming  and  for  this  par¬ 
ticular  case,  the  ASR  of  the  membrane  should  be  less  than 
7  £2  cm2.  Most  of  the  reforming  takes  place  within  the  first 
one-third  of  the  reformer  channel.  Also  shown  in  Fig.  2  is 
the  reformate  concentration  profile  in  the  anode  channel  for 
external  reforming.  This  is  strictly  linear  because  the  con¬ 
stant  cell  current  removes  reformate  at  a  constant  rate.  The 
corresponding  anode  profile  for  membrane  reforming  shows 
a  slight  curvature  caused  by  the  membrane  increasing  the 
reformate  concentration  in  the  anodes  of  the  low  index  cells. 

Fig.  3  shows  the  oxygen  activity  profiles  in  the  reformer 
and  anode  channels.  It  is  the  difference  between  them  that 
drives  permeation  through  the  reformer  membrane.  The  dif¬ 
ference  is  largest  for  the  lower  index  cells,  which  is  consis¬ 
tent  with  most  of  the  reforming  being  performed  in  these 
cells.  The  oxygen  activity  profiles  in  the  anode  are  very  sim¬ 
ilar  for  both  the  membrane  and  external  reforming  cases,  as 
is  expected  from  the  similarity  in  the  reformate  profiles  in 
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Fig.  2.  Profiles  of  the  fraction  of  reformate  in  the  gas  stream  in  the  reformer  channel  and  the  anode  channel  calculated  from  the  membrane  reformer 
model  at  the  critical  membrane  ASR  of  7  £2  cm2  (cell  current  density:  0.25  A  cm-2,  temperature:  900  °C,  fuel  utilisation:  85%).  Also  shown  is  the  profile 
in  the  anode  channel  for  external  steam  reforming  under  the  same  conditions. 


Fig.  2.  The  cell  voltage  profiles  are  also  shown  in  Fig.  3 
and  there  is  no  detectable  difference  between  the  membrane 
and  external  reforming  cases.  The  driving  force  for  trans¬ 
port  across  the  reformer  membrane  can  be  expressed  as  an 


electro-motive  force  given  by  applying  the  Nernst  equation 
to  the  ratio  of  oxygen  partial  pressures  on  the  two  sides 
of  the  membrane.  This  is  also  shown  as  a  function  of  cell 
number  in  Fig.  3. 
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Fig.  3.  Profiles  of  oxygen  activity  in  the  reformer  channel  and  the  anode  channel  calculated  from  the  membrane  reformer  model  at  the  critical  membrane 
ASR  of  7  £2  cm2  (cell  current  density:  0.25  A  cm-2,  temperature:  900  °C,  fuel  utilisation:  85%).  Also  shown  is  the  profile  in  the  anode  channel  for 
external  steam  reforming  under  the  same  conditions,  the  cell  voltage  profile  for  both  conditions  (these  are  indistinguishable)  and  the  Nernst  emf  across 
the  membrane. 
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Table  1 


Summary  of  calculations  comparing  a  SOFC  with  membrane  reforming 
with  one  having  external  steam  reforming 


Parameter 

External 

reforming 

Membrane 

reforming 

Membrane  ASR  (£2  cm2) 

N/A 

7 

Fuel  utilisation  (%) 

84.2 

84.9 

Electrical  efficiency  (%  LHV) 

48.2 

62.4 

Total  voltage  (V) 

16.0 

15.89 

Fuel  input  velocity  (ms-1) 

0.119 

0.031 

Exhaust  output  velocity  (ms-1) 

0.119 

0.0892 

Operating  temperature:  900  °C,  current  density:  0.25  A  cm  2,  SOFC  ASR: 
0.5  Q  cm2. 


The  electrical  efficiency  was  calculated  for  the  SOFC  with 
membrane  reformer,  the  output  electrical  power  divided  by 
the  room  temperature  LHV  of  the  input  methane  (A #298  = 
—802,  286  J mol-1  [11]).  This  was  also  done  for  the  exter¬ 
nal  reforming  base  case.  The  results  are  shown  in  Table  1. 
The  membrane  reforming  increases  the  efficiency  for  elec¬ 
trical  power  generation  from  48.2  to  62.4%.  This  is  entirely 
due  to  the  fact  that  the  internal  membrane  reforming  essen¬ 
tially  converts  the  SOFC  to  direct  methane  operation.  The 
lower  efficiency  of  external  steam  reforming  is  due  to  the 
enthalpy  requirement  of  the  steam  reforming  reaction.  If  all 
this  enthalpy  could  be  extracted  from  the  waste  heat  of  the 
SOFC  then  the  efficiency  for  steam  reforming  would  in¬ 
crease  to  that  of  the  membrane  reforming  case.  Since  inter¬ 
nal  steam  reforming  should  be  able  to  achieve  this,  there  is 
no  efficiency  advantage  of  membrane  reforming  over  inter¬ 
nal  steam  reforming,  but  membrane  reforming  removes  the 
other  problems  associated  with  on-anode  steam  reforming. 


4.  Potential  membrane  material 

The  family  Lai-^Sr^Fei-yCr^Os  (abbreviated  as  LSFCr) 
contains  potentially  suitable  MIEC  compositions  whose  sta¬ 
bility,  mass  transport  and  interfacial  reaction  rates  can  be  ad¬ 
justed  by  changing  v  and  y  [12,13].  Increasing  v  or  decreas¬ 
ing  y  increases  mass  transport,  but  reduces  stability  (both 
chemical  and  mechanical).  The  best  compromise  was  found 
to  be  with  y  =  0.2  and  v  in  the  range  0.4-0. 6.  Oxygen  tracer 
diffusion  coefficients  and  isotopic  surface  exchange  con¬ 
stants  measured  under  reducing  conditions  similar  to  those 
predicted  for  membrane  reforming  are  summarised  in  Fig.  4 
[14].  At  900  °C  and  under  reducing  conditions  the  oxygen 
tracer  diffusion  coefficient  is  approximately  10_7cm2s_1 
and  is  relatively  insensitive  to  oxygen  activity  [15].  Thus 
the  integration  in  Eq.  (2)  is  not  necessary.  From  Eq.  (3)  for 
a  membrane  thickness  of  0.5  mm  this  tracer  diffusion  co¬ 
efficient  is  equivalent  to  a  contribution  of  1 .6  £2  cm2  to  the 
ASR.  Since  the  target  value  is  <7  £2  cm2,  diffusion  in  these 
materials  is  sufficiently  fast  for  this  application. 

Under  the  same  conditions  the  surface  isotopic  exchange 
rate  constant  with  H2O/H2  (&*)  is  in  the  range  3  x  10~7 


10-5  7 


■  •  X  =  0.2 

■  T  x  =  0.4 
■  x  =  0.6 

10‘8  -I - 1 - 1 - 1 - 1 - 1 - 

0.75  0.80  0.85  0.90  0.95  1.00  1.05 

103/T  (K"1) 

Fig.  4.  D*  and  k  measured  for  LSFCr  compositions  in  a  gas  composition 
of  0.05  atm  H2,  0.15  atm  N2,  0.1  atm  H2O. 

to  10  x  10-7  cms-1  at  a  water  vapour  pressure  of  0.1  atm. 
The  reaction  rate  constant  is  related  to  the  tracer  exchange 
constant  by  [16] 

k  =  k*  (2 n  +  A  (20) 

where  n  is  the  order  of  reaction  for  the  gaseous  species,  y 
is  a  thermodynamic  factor 

1  d(ln  Po?) 

y  = - 1 - 'C1L  (21) 

r  2  d(lnco) 

where  co  is  the  concentration  of  oxygen  in  the  solid  phase. 
The  dependence  of  k*  on  water  vapour  pressure  is  not  known, 
but  could  be  linear  i.e.  n  =  1  (it  is  assumed  to  be  independent 
in  the  model).  Similarly,  y  is  not  known  but  is  likely  to 
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be  ^1  in  the  region  where  D*  is  independent  of  oxygen 
activity.  Thus  Eq.  (20)  reduces  to  k  =  2 k*.  The  water  vapour 
pressure  varies  along  the  anode  channel,  but  that  used  in 
the  measurements  (0.1  atm)  is  a  typical  value.  Hence  from 
Eq.  (5)  the  contribution  to  the  membrane  ASR  from  the 
cathodic  reaction  with  water  is  in  the  range  1.6-5. 3  £2  cm2 
at  900  °C. 

At  the  surface  of  the  membrane  with  the  reformer  channel 
the  anodic  reaction  parameter  ka  is  the  reaction  rate  constant 
for  oxidation  of  methane  by  LSFCr.  The  turnover  frequency, 
Nt,  for  this  reaction  has  been  measured  to  be  approximately 
three  molecules  of  CH4  per  site  per  second  at  5%  methane 
and  900  °C  [17].  The  experiments  also  showed  that  these 
materials  are  extremely  resistant  to  carbon  deposition  even 
when  exposed  directly  to  dry  methane.  The  rate  constant,  as 
defined  by  Eq.  (3),  is  related  to  the  turnover  frequency  by 


NtV\ mo  RT 
NaM  A/xo 


(22) 


ing  steam  and  anode  cooling  encountered  with  conventional 
on-anode  steam  reforming.  The  simulations  indicate  an  in¬ 
crease  in  electrical  efficiency  of  approximately  14  percent¬ 
age  points  with  respect  to  external  steam  reforming  at  the 
same  fuel  utilisation.  For  a  fuel  cell  ASR  of  0.5  £2  cm2,  the 
membrane  should  have  an  ASR  of  less  than  7  £2  cm2. 

Experimental  measurements  of  oxygen  diffusion  and  sur¬ 
face  reaction  show  that  materials  in  the  family  Lai_xSrx 
Fei-yCvyOi  have  sufficiently  fast  oxygen  diffusion  under 
reducing  conditions  to  meet  the  target  ASR.  However,  the 
surface  reaction  rates  (especially  oxidation  of  methane  at 
the  reformer  surface  of  the  membrane)  are  probably  too  low 
by  about  a  factor  of  2.  The  rates  could  be  increased  by  us¬ 
ing  porous  coatings  of  the  same  material  on  the  membrane 
surfaces  in  order  to  meet  the  target. 
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where  As  is  the  area  of  a  site  measured  by  nitrogen  adsorp¬ 
tion  (1.62  x  10-19m2).  The  difference  in  oxygen  chemi¬ 
cal  potential  driving  the  reaction  in  the  measurements  of  Nr 
was  not  known,  but  the  driving  force  was  large  and  proba¬ 
bly  beyond  the  limit  of  linear  behaviour  on  which  Eq.  (3)  is 
based.  In  order  to  make  an  approximate  estimate,  we  there¬ 
fore  assume  the  factor  RT/ A/xo  to  be  of  order  unity.  The 
turnover  frequency  is  proportional  to  the  methane  concen¬ 
tration  (although  assumed  independent  in  the  model)  and  in 
the  reformer  channel  this  is  approximately  0.5  atm  which  is 
10  x  that  used  in  the  measurement  of  Nr.  Hence  Nr  for  the 
anodic  membrane  reaction  is  approximately  30s_1  per  site 
and  the  corresponding  value  of  ka  for  methane  oxidation  is 
thus  3.7  x  10-7  cms-1.  From  Eq.  (5)  this  gives  a  contribu¬ 
tion  to  the  membrane  ASR  of  8.6  £2  cm2. 

The  total  membrane  ASR  is  thus  approximately  in  the 
range  12-16  £2  cm2  which  is  about  double  the  target  value 
from  the  membrane  reformer  model.  The  experimental  data 
indicate  that  the  surface  reactions  dominate  the  equivalent 
resistance  of  the  membrane  and  that  the  largest  equivalent 
resistance  is  probably  the  reaction  with  methane.  Thus,  the 
membrane  surfaces  are  likely  to  require  catalytic  coatings  to 
increase  the  interfacial  reaction  rates.  However,  the  required 
increases  are  modest  (less  than  an  order  of  magnitude)  and 
this  could  probably  be  achieved  by  using  a  porous  coating 
of  the  same  LSFCr  membrane  material. 


5.  Conclusions 

A  simplified  reaction  model  has  shown  that  in- stack  mem¬ 
brane  reforming  in  a  SOFC  is  conceptually  feasible.  Mem¬ 
brane  reforming  has  the  same  efficiency  as  operating  the 
SOFC  directly  on  dry  methane,  but  with  much  less  risk  of 
carbon  deposition.  It  also  avoids  the  complications  of  rais- 


The  authors  are  grateful  to  EPSRC  and  the  Link  Ap¬ 
plied  Catalysis  Program  for  financial  support  under  grant 
GR/M05928.  T.  Ramos  is  also  grateful  to  FCT-Programa 
Praxis  XXI  (Portugal)  for  the  award  of  a  study  scholarship. 

References 

[1]  H.  Kim,  C.  Lu,  W.L.  Worrell,  J.M.  Vohs,  R.J.  Gorte,  J.  Electrochem. 
Soc.  149  (2002)  A247. 

[2]  T.  Hibino,  S.Q.  Wang,  S.  Kakimoto,  M.  Sano,  Solid  State  Ionics  127 
(2000)  89. 

[3]  P.  Aguiar,  E.  Ramirez-Cabrera,  N.  Laosiripojana,  A.  Atkinson,  L.S. 
Kershenbaum,  D.  Chadwick,  Stud.  Surf.  Sci.  Catal.  145  (2003)  387. 

[4]  P.  Vernoux,  M.  Guillodo,  J.  Fouletier,  A.  Hammou,  Solid  State  Ionics 
135  (2000)  425. 

[5]  B.C.H.  Steele,  in:  H.L.  Tuller  (Ed.),  Oxygen  Ion  and  Mixed  Conduc¬ 
tors  and  Their  Technological  Applications,  Kluwer  Academic  Pub¬ 
lishers,  Dordrecht,  2000,  p.  323. 

[6]  U.  Balachandran,  J.T.  Dusek,  R.L.  Mieville,  R.B.  Poeppel,  M.S. 
Kleefisch,  S.  Pei,  T.P.  Kobylinski,  C.A.  Udovich,  A.C.  Bose,  Appl. 
Catal.  A:  General  133  (1995)  19. 

[7]  J.E.  tenElshof,  M.H.R.  Lankhorst,  H.J.M.  Bouwmeester,  J.  Elec¬ 
trochem.  Soc.  144  (1997)  1060. 

[8]  S.  Kim,  Y.L.  Yang,  A.J.  Jacobson,  B.  Abeles,  Solid  State  Ionics  121 
(1999)  31. 

[9]  E.F.  Sverdrup,  C.J.  Warde,  R.L.  Eback,  Energy  Conv.  13  (1973)  129. 

[10]  F.J.  Gardner,  M.J.  Day,  N.P  Brandon,  M.N.  Pashley,  M.  Cassidy,  J. 
Power  Sour.  86  (2000)  122. 

[11]  HSC  Chemistry  v4,  Outokumpu,  Pori,  Finland. 

[12]  W.T.  Stephens,  T.J.  Mazanec,  H.U.  Anderson,  Solid  State  Ionics  129 
(2000)  271. 

[13]  A.  Atkinson,  T.  Ramos,  Solid  State  Ionics  129  (2000)  259. 

[14]  T.  Ramos,  A.  Atkinson,  in:  T.A.  Ramanarayanan,  W.L.  Worrell, 
M.  Mogensen  (Eds.),  Ionic  and  Mixed  Conducting  Ceramics  IV, 
Electrochemistry  Society,  2001. 

[15]  T.  Ramos,  Ph.D.  Thesis,  Imperial  College,  London,  2002. 

[16]  S.  Kim,  Y.L.  Yang,  A.J.  Jacobson,  B.  Abeles,  Solid  State  Ionics  121 
(1999)  31. 

[17]  E.  Ramirez-Cabrera,  Ph.D.  Thesis,  Imperial  College,  London,  2001. 


